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A new Monte Carlo study of evolution 



We solve CCFM evolution equation numerically using the CohRad program based 
on Monte Carlo methods. We discuss the effects of removing soft emissions and 
non-Sudakov form factor by comparing the obtained distributions as functions of 
accumulated transverse momenta or fractions of proton's longitudinal momenta. We 
also compare the solution of the CCFM with the DGLAP equation in the gluonic 
channel. Finally, we analyze the infra-red behaviour of solutions using the so-called 
diffusion plots. 
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Abstract 



1 Introduction 



The Large Hadron Collider opened up the possibility to scan parton densities over a wide 
domain of partons kinematics. This allows for detailed studies of various theoretically in- 
teresting and phenomenologically relevant dynamical effects taking place during partons 
evolution, such as: coherence yj, saturation |2| or both (3]-[5]. In the present work we are 
in particular interested in the coherence effects in the initial state gluon cascade which 
is modelled by CCFM [l||6) evolution equation. These effects are taken into account by 
summing up dominant contribution in angular ordered regions of phase space. There- 
fore, in comparison with DGLAP |7| equation, CCFM includes some of the interference 
effects that are subleading from the point of view of approximation to leading-order in 
the ordering in the hard scales In where Qt is the transverse momentum of a t 

channel gluon. Because of this type of ordering, the CCFM equation is applicable also in 
the domain of low x and can be viewed as a bridge between low x and large x physics. 
The CCFM equation, due to the coherence effects included, provides parton distribution 
function (PDF) not only at a given fraction of proton longitudinal momentum x and 
transverse momentum accumulated in the gluonic ladder, but accounts also for additional 
argument p, related to the maximal angle of gluon emissions. Hence, it allows for match- 
ing a PDF with a hard process matrix element given the scale of the last emission. The 
last feature makes it particularly interesting for extextension of the BFKL fs] approach. 
The classical CCFM equation is linear and predicts unlimited growth of parton densities 
at small x. It can, however, be extended in order to account for saturation by extending it 
by a non-linear term |3|[9| or impose absorptive boundary conditions 10 IT]. The CCFM 
equation has been already studied theoretically in some limiting cases i.e. in the low x 
limit jT0,12-T4] and within Monte Carlo formulation 15 -18). It has been also used in 



phenomenological applications 22 . However, the open issues concerning the CCFM 
are numerous to mention here just the proper form of the initial conditions, the details of 
violation of unitarity, the role of the soft emissions in the low x limit. Also the efficient 
algorithm for solving it together with evolution of quark densities is still an open problem. 
In the present study we reinvestigate systematically aspects of the CCFM equation using 
the Monte Carlo Markov chain approach focusing on the role played by the non-Sudakov 
form factor and the soft parts of the splitting function. The understanding of the effects 
coming from both parts of the splitting functions is important for investigations of uni- 
tarity violation effects in evolution of partons |5) since the nonlinearities affect the soft 
emissions. Such investigation is possible since the Monte Carlo numerical integration we 
use [23j allows for easy handling of singular integrals and therefore for taking into account 
effects from the full splitting function i.e. soft and hard emissions and corresponding form 
factors. The second reason linked to the use of Markov chain Monte Carlo is to provide 
a new scheme for performing numerically efficient parton shower based on a forward evo- 
lution and extend it in the future to account for a nonlinear term allowing for saturation 
of gluons as well as for quarks. 

The paper is organized as follows. In section 2 we introduce an iterative formulation of the 
CCFM equation suitable for the Markov Chain implementation. In section 3 we present 
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the Monte Carlo algorithm we apply to solve the CCFM equation using the CohRad 
program. In the section 3 we analyze properties of the CCFM equation i.e. the effect of 
neglecting the Sudakov form factor and soft real emissions, and we compare the CCFM to 
angular ordered DGLAP cascade. Finally we perform the analysis of population of differ- 
ent regions of phase space focusing on the diffusion aspects of the considered equations. 

2 CCFM evolution equation 




Figure 1: Notation used in the text. 

We use notation as in Fig. [T] pi and qi denote four-momenta of real and virtual 
gluons, respectively, and Xj+i = Zj+iXj are fractions of longitudinal momenta of the gluon 
initiating the cascades. If 2; ~ 0, a large momentum fraction has been carried out by a 
real emitted gluon (hard emission), while 2 ~ 1 corresponds to a soft emission. 2-vectors 
of transverse momenta of the emitted gluons are denoted by p^rp and transverse momenta 
accumulated on the emitting line by qiT = |qoT — Yl]=i Pjt|- It is also convenient to use 
rescaled transverse momenta pi = with their modulus being related to angles of 

emissions: pi = |pj| = ~ -Ej^j- 

The CCFM equation imposes angular ordering in the real emissions, that can be 

expressed either using angles: pi^i > Zipi or rapidities of the emitted gluons: In — = 

ViT 

rji < Vi+i- 

In the following we present solution of the CCFM equation for the unintegrated gluon 
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density in the iterative form: 



A{x,qT,p) = A{xo,qoj.,Po)^s{p,Po)S{x - Xo)S{qT - qoj 
+ ^ dxo (fqoTA{xo,qoT,Po) 
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convenient for a Markov chain Monte Carlo implementation. The Q{pi — Zi_iPi_i) func- 
tions impose angular ordering of emissions. The scale p related to rapidity position of 
the hard process will be defined more precisely in the following. The variable po = 1 
GeV plays the role of the minimal scale and the infrared cutoff on transverse momenta, 
PiT > Po- The initial transverse momentum of a gluon coming from a proton is denoted 
by goT- In the above 

p / ~ N (^sNc f 1 ANs{p,z,qT)\ 

Pgg{z,P, <1t) = — [—, + ) (2) 

is the CCFM splitting function. Its form is similar to the DGLAP Q splitting: 

nr^''w = ^(r^ + -2) (3) 

yet does not include terms non-singular in z — )■ 0. The non-Sudakov form factor A7V5'(p, z, g^) 
reads: 

A.,(p,.,,,) = exp|-— -y^^^_^-^|=exp|-— In-ln — |. (4) 

The particular form of the non-Sudakov form factor we study here and which allows 
the form factor to be larger than unity was motivated by the investigations in [ll] . This 
formulation which we wanted to reproduce in Monte Carlo neglects kinematical constraint 
effects and leads to fast growth of gluon density towards small values of qx- Studies taking 
into account kinematical constraint effects we postpone for future investigations. In both 
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Arts and as well as in the splitting function we kept constant, = 0.2, for simplicity 
at this stage. In the future we plan to use coupling constant as suggested by NLO BFKL 
results |24j. This form factor enters the CCFM equation from resummation of virtual 
emissions that are harder than either of the emitting lines. It regulates z — )■ singularity 
of the splitting function. A detailed discussion on the physical interpretation of the region 
of integration can be found for instance in ref. [lOj . 
The Sudakov form factor is given by 



A,(g,g')=exp ^ . (5) 




It can be interpreted as resummed soft virtual emissions on the emitter line. With po 
being the minimal allowed transverse momentum|^ the soft (IR) cutoff in CCFM e{pi) = 
Po/fexo) is evolution scale dependent. At the beginning of the evolution {pi ~ po) only 
hard, Zi — )■ 0, emissions are allowed. In DGLAP this cutoff is constant, e{pi) = e. 

Finally let us also mention that the distribution A in eq. ([T| is related to the gluon 
density g{x) through the relation 

xg{x,Q')= [ ^xA{x,qT,Q')Q{Q'-q'T). 



3 Monte Carlo algorithm 

We implemented CCFM evolution equation ([T| in CohRad as a Markov chain in (r]i,Xi). 
Due to the similarities between the structures of CCFM and DGLAP, we applied impor- 
tance sampling and generated emissions from DGLAP distributions that were reweighted 
to obtain CCFM. This is a standard Monte Carlo technique, that improves efficiency of 
generation. In the DGLAP evolution the Sudakov form factor has been used as a prob- 
ability distribution in generating each emissions and as a correcting weight for the last 
emission. The details of the algorithm can be found in ref. |25| . 

The emitted particles' four-momenta are constructed from rji and Xi using the follow- 
ing: 

Pj^ = exp{-r]j)^{xj^i - Xj), p+ = Vs{xj_i - Xj), pj = (6) 

pj 

It can be checked that rjj^i > rjj hence angular ordering of emissions is imposed by 
construction. The maximal allowed rapidity in the hard process frame is equal to 0. Its 
smallest value is determined by the infrared cutoff on real emissions. We set it on the 
minimal transverse momentum pQ = p™" = IGeV. The rapidity of the first emitted gluon 
must therefore obey 

Po < e-'^^y^xo, (7) 
^ We set Po =lGeV, see below for more discussion. 
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which leads to Xo-dependent minimum rapidity 




(8) 



The sphtting function used in the Monte Carlo is 
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(9) 



The non-Sudakov form factor A/vr5(p, z, q) is included in the form of a correcting weight 
for a "DGLAP-wise" generated distributions. It multiplies both 1/z and 1/(1 — z) poles. 

4 Results 

We present solutions of CCFM evolution equation for unintegrated gluon distributions 
xA{x,qT,p)- They are functions of three parameters: fractions of proton's longitudinal 
momentum x, transverse momenta Qt accumulated on the emitting line and the hard scale 
p. The role of hard scale p in eq. ([T| is to provide an upper limit for pi. 

Let us first compare in more detail the implementations of the CCFM and DGLAP 
evolution equations, before moving to analyzing numerical results. The differences be- 
tween DGLAP and CCFM lie first of all in the absence of non-Sudakov formfactor in 
DGLAP. Both evolution equations include Sudakov form factor, yet differ in the defini- 
tion of infrared cutoff e, see discussion after equation ([s]). The scale-dependent IR cutoff 
in CCFM suppresses the evolution into low values of x for small scales, whereas DGLAP 
allows for emissions in that region. Another important point is the definition of the hard 
scale p in case of CCFM and DGLAP. In case of CCFM, following ref. [Is] , it is set to fixed 
(x-independent) angle, p = Q. In terms of maximal rapidity, however, the above leads to 
ln(Q/x), see also ref. fTs]. The variable x = Xn refers to the gluon entering the hard pro- 
cess. Concerning the hard scale pi, in our implementation of DGLAP (following ref. |25j) 
we have the hard scale set as the maximal rapidity ln(Q) and it is x-independent. That 
translates into maximal angles, that are x-dependent p{x) = xQ. Due to the dependence 
of p on the final x, the eq. ([T]) gets slightly modified: 



The above is easily implementable in the Markovian type of MC. It is this effect that 
causes CCFM distribution to rapidly rise at small values of x. 

It is important that for DGLAP the Monte Carlo program enables us to extract not 
only X but also qt and treat the solution of the DGLAP evolution equation in the following 
numerical analysis as a function of three arguments x, qx and p. The initial distribution 
will be typically XQAQ{xo,qTQ,p) = ^/Qtq both, CCFM and DGLAP 




(10) 



Pi<p 



Pi <p{Xi-l ) 
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Figure 2: Comparison of distributions in x (integrated over transverse momenta) for 
different hard scales: Q =100GeV (left) and lOOOGeV (right). 

In Fig. [2] the two solutions are shown as functions of x (integrated over transverse 
momenta) at different hard scales Q. In both cases the starting condition is uniform in 
X in order to observe the effect of evolution itself on the distributions. It is visible that 
CCFM equation generates a much sharper rise in the small values of x. This is caused 
by the hard scale p in CCFM that allows the emitting more gluons, due to larger (by 
ln{xQ/x)) allowed rapidity range than in DGLAP. As a result, more emissions occur, and 
each of them takes away more longitudinal momentum than in the case of DGLAP. 




loglO(qT) loglO(qT) 

Figure 3: Comparison of solutions obtained in CCFM (solid lines) and DGLAP evolution 
(dashed lines). Two plots were made for different maximal angle Q = 100, lOOOGeV. Initial 
condition: 

In Fig. [3] we compare CCFM and DGLAP distributions of accumulated transverse 
momenta Qt- Different curves correspond to different values of final x, as indicated in the 
plots. The relative differences in normalization for different regions in x come from the 
rescaling factors . . — r^. In the region below qr = IGeV both distributions reflect 

X- width 01 X bin ° 

the initial condition. The kink at qt = IGeV is related to the infrared cutoff on the 
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minimal transverse momenta of emitted particles. Above IGeV the shapes of distributions 
are generated by evolutions. It can be seen that at large values of x both distributions 
are equivalent. At small values of x CCFM distribution falls down with increasing qt 
steeper than DGLAP. This suppression can be attributed to the non-Sudakov form factor 
which is smaller than 1 for sufficiently large-angle emissions: pi > gj^/y^. In order to 
show that the shape of distributions is influenced by evolution only in the region of large 
transverse momenta, the simulations were repeated with initial conditions exp(— g^). The 
results are shown in Fig. |4j 
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Figure 4: Left: CCFM evolution with initial conditions: (solid lines) and exp{—q^) 

(dashed lines). Right: The same plot for DGLAP solution. The hard scale in both plots 
is g = lOOOGeV. 




Figure 5: Solution of complete CCFM equation (solid lines) and equation with soft emis- 
sions removed. The hard scale is: Q = lOOGeV (left), lOOOGeV (right). Initial condition 

The CCFM equation is often solved in a simplified form, with the Sudakov form factor 
and 1/(1 — 2;) pole in the splitting function removed. This approximation, described for 



7 



instance in 11 , 12 , is valid in the small x region and for inclusive observables. Fig. [Sjshows 
comparison between distributions obtained from a complete and simplified equations for 
two maximal angles and two values of x. This approximation appears to work well, giving 
similar slope of the distributions in the region of large transverse momenta for small x 
(dark blue lines). The similarity of the distributions at very small values of Qt reflect 
only the common initial condition. For large values of x the complete and simplified 
equations diverge, as expected. Neglection of the soft emissions i.e. Sudakov form factor 
and 1/(1 — 2;) pole manifests itself particularly strongly at low qx since we observe a sharp 
cliff on Fig. 5 in small values of qt ~ 10^^. In the Monte Carlo simulation this reflects 
itself as small MC weights given to soft events. 

Another interesting aspect of comparing DGLAP and CCFM equations lies in observ- 
ing the way their solutions populate the region corresponding to infrared momenta. In 
Figs. (Gllslwe show two-dimensional distributions 'i!{x, qT,p) obtained from solutions of the 



evolution equation A{x,qT,p) 

,T,/ qTxA{x,qT,p) , . 

^{x,qT,p) = 77 r, 11 

V 'yj'/^v qjp''''{x)xA{x,qrp''='{x),p)' ^ ' 

where qjp°-^{x) is the value of transverse momentum, for which qTxA{x, q,TP) is maximal. 
In the following, we restrict transverse momenta to the values larger than 1 in order to 
study the effects of evolution only. The form of \l/(a:, qT,p) is motivated by the asymptotic 
behaviour of (linear) BFKL equation, see ref. (26) for details. 




Figure 6: The diffusion ratio ^ of eq. (11) for CCFM, Q = 100 (left), 1000 (right). 



The diffusion plots in Figs. iGjjSj are made according to convention of ref. |26| for all 
three considered evolutions and hard scales Q = 100, 1000 (left, right plots, respectively). 
Fig.|6]shows distribution obtained from the complete CCFM equation, Fig.[7]from DGLAP 
evolution, and finally Fig. [8] from CCFM without soft emissions. On the plots we only 
show the distribution in the phase space region > IGeV. Comparing Fig. 6 and Fig. 
8 we see that at low x region the solution of the CCFM is well approximated by just hard 
emissions. By comparing Fig. 7 to Figs. 6 and 8 we see that the presence of non-Sudakov 
leads to effectively more diffusive distribution in what is consistent with BFKL like 
properties of CCFM. 
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Figure 7: The diffusion ratio \1/ of eq. (11) for DGLAP. 





Figure 8: Tlie diffusion ratio \1/ of eq. (11) for CCFM witliout soft emissions 



5 Conclusions 



In tliis paper we studied some of tlie properties of the CCFM evolution equation formu- 
lated as a Markov chain parton shower model. This approach is particularly interesting 
since it can be generalized to a simulation for a nonlinear evolution equation, and will 
be presented elsewhere. The purpose of the paper was to demonstrate capabilities of 
this method as well as to check some of the properties of CCFM which are difficult to 
quantify by other methods e.g. neglecting contributions on inclusive level of Sudakov 
form factor and soft emissions at small x. In the standard non Monte Carlo calculations 
it is a formidable task to perform a numerical calculation with both 1/z and 1/(1 — z) 
poles included. Our calculation supports the claim that at low x one can safely neglect 
contribution of both Sudakov form factor and soft emissions i.e. 1/(1 — z) pole. Another 
subject we studied concerned the comparison of the CCFM dynamics with DGLAP type 
of resummation. We observed that the CCFM equation due to angular ordering and par- 
ticular condition for the maximal angle of emitted partons generates much faster growth 
of gluon density than DGLAP at small x. Finally we addressed the question how the 
partons generated according to CCFM and DGLAP equations populate the {x^qt) plane. 
The last study will be useful for future investigations of regions of phase space where 
the unitarity is violated bringing the need for suppression of growth of gluons by their 
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recombination. 
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